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We show that irradiation of a voltage-biased superconducting quantum point con- 
tact at frequencies of the order of the gap energy can remove the suppression of 
subgap dc transport through Andreev levels. Quantum interference among reso- 
nant scattering events involving photon absorption is furthermore shown to make 
microwave spectroscopy of the Andreev levels feasible. We also discuss how the 
same interference effect can be applied for detecting weak electromagnetic signals 
up to the gap frequency, and how it is affected by dephasing and relaxation. 
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1. Introduction 

It is well known that the current through a voltage-biased superconducting quantum point contact 
(SQPC) is carried by localized states. These states, called Andreev states, are confined to the normal 
region of the contact. The energy of the states — the Andreev levels — exist in pairs, (one above and 
one under the Fermi level), and lie within the energy gap of the superconductor, with positions which 
depend on the change </> in the phase of the superconductors across the junction. The applied bias 
affects this phase difference through the Josephson relation, (j> — 2eV/h. With a constant applied 
bias V much smaller than the gap energy A, <j) will increase linearly in time, and the Andreev levels 
will move adiabatically within the gap. This motion is a periodic oscillation in 0, indicating that 
no energy is transfered to the SQPC and a pure ac current will flow through the contact. This is 
actually the ac Josephson effect. 

We wish to study this system in a non-equilibrium situation, one way to accomplish this is by 
introducing microwave radiation with a frequency to « A, which will couple the Andreev levels 
to each other. The radiation will represent a non-adiabatic perturbation of the SQPC system. 
However, if the amplitude of the electromagnetic field is sufficiently small, the field will not affect 
the adiabatic dynamics of the system much unless the condition for resonant optical interlevel 
transitions is fulfilled. Such resonances will only occur at certain moments determined by the time 
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Fig. 1. Time evolution of Andreev levels (full lines) in the energy gap of a gated voltage-biased, 
single-mode SQPC (see inset). A weak microwave field induces resonant transitions (wavy lines) 
between the levels at points A and B and the level above the Fermi energy becomes partly occu- 
pied to an extent determined by interference between the two transition amplitudes. Non-adiabatic 
interactions release the energy of quasiparticles in the (partly) occupied Andreev level into the con- 
tinuum at point C, where the Andreev states and the continuum merge into each other (represented 
by dashed arrows, see text) and the initial conditions for the Andreev level populations are reset 
(filled and empty circles). 



evolution of the Andreev level spacing. The resonances will provide a mechanism for energy transfer 
to the system to be nonzero when averaged over time and hence for a finite dc current through the 
junction. The rate of energy transfer is in an essential way determined by the interference between 
different scattering events ||, and therefore it is not surprising that oscillatory features appear 
in the (dc) current-voltage characteristics of an irradiated SQPC. Dephasing and relaxation will 
affect the interference pattern and may even conspire to produce a dc current flowing in the reverse 
direction with respect to the applied voltage bias. The mechanism behind this negative resistance 
is very similar to the one responsible for the "somersault effect" discussed by Gorelik et al. j|. 

A peculiar feature of the Andreev bound states, in comparison with normal bound states, is that 
they can carry current. This is why microwave-induced transitions between Andreev levels can be 
detected by means of transport measurements. In fact, it is possible to do Andreev energy-level 
spectroscopy in the sense that the Andreev level spectrum at least in principle can be reconstructed 
from a measurement of the microwave-induced subgap current. Such microwave spectroscopy of 
Andreev states in mesoscopic superconductors is the topic addressed in this work. 



2. Theoretical framework 

For an unbiased, mesoscopic SQPC — which by construction has a normal region L which is much 
shorter than the coherence length £o — the Andreev spectrum of each transport mode has the form 

E ± ((f>) = ±E{<f) = ±A[1 - D sin 2 (0/2)] Va. (l) 

Here D is the transparency of the mode and the energy is measured from the Fermi energy Q J5|. 
With a small bias voltage applied, the levels move along the adiabatic trajectories E ± (t) = ±E((f> + 
2eVt/K) in energy-time space, oscillating with a period of T p = fnr/eV, as shown in Fig. I. In 
equilibrium the lower Andreev level will obviously be occupied, while the upper level will be empty. 
In the discussion below we consider a single-mode SQPC, although we note that the theory also 
applies for the case of a single dominant mode in a multi-mode junction. The transparency D of the 
mode is taken to be arbitrary but energy-independent, < D < 1. 

A high frequency electromagnetic field is applied to the gate situated near the contact, see inset 
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in Fig. 1. The time-dependent electric field induced by the gate is concentrated within the non- 
superconducting region of the junction, and hence the charge carriers will couple to the electromag- 
netic field only there. When the criterion Txj) = 2eV <C 2E 2 (t)/A for adiabaticity is obeyed, the rate 
of interlevel transitions is exponentially small keeping the level populations constant in time || ^|. 
The presence of a weak electromagnetic field [on the scale of E(t)] does not affect the adiabatic level 
trajectories except for short times close to the resonances at t = t^ b, when E(t^B) — fuo/2. Here 
the dynamics of the system is strongly non-adiabatic with a resonant coupling which effectively 
mixes the adiabatic levels. This is an analog of the well known Landau-Zener transition, which de- 
scribes interlevel scattering as a resonance point is passed. In our case these transitions give rise to 
a splitting of the quasiparticle trajectory at the points A, B into two paths; A1A2B2 and A\B\B2, 
forming a loop in (E,t) space (see Fig. 1). 

The resonant scattering opens a channel for energy absorption by the system; a populated upper 
level when approaching the edge of the energy gap (at point C in Fig. 1) creates real excitations in 
the continuum spectrum, which carry away the accumulated energy from the contact. As a result, 
the net rate of energy transfer to the system is finite; it consists of energy absorbed both from the 
electromagnetic field and from the voltage source. The confluence of the two adiabatic trajectories 
at B2 gives rise to a strong interference pattern in the probability for real excitations at the band 
edge, point C. The interference effect is controlled by the difference of the phases acquired by the 
system during propagation along the paths A1A2B2 and A\B\B2- 

2.1. Bogoliubov-de Gennes equation 

In order to describe the time evolution of the Andreev states in more detail we use the time- 
dependent Bogoliubov-de Gennes, (BdG), equation B for the quasiclassical envelopes u±(x,t) of 
the two-component wave function ^(x, t) = u+(x, t)e lkFX + u_ (x, t)e pX , 

iKdu/dt=[H + o- z V g (x 1 t)]u. (2) 

In this equation u = (u + ,u_) is a four-component vector, while Ho is the Hamiltonian of the 
electrons in the electrodes of the point contact, 

H = — ihvpa z T z d/dx + A{cos(<p(t)/2)a x + i sin((p(t) /2)sgaxa y ], (3) 

where <7j and n denote Pauli matrices in electron-hole space and in u± space respectively. For a 
mesoscopic junction the gap function A(r) need not be calculated self-consistently, which is why in 
Eq. (g) it is assumed to have the constant value A exp(i(/)/2) in the superconducting reservoir on one 
side of the SQPC and the different constant value Aexp(— i<j)/2) in the reservoir on the other side. 
The gate potential V g {x,t) = V UJ (x) cosujt in Eq. (||) oscillates rapidly in time and the amplitude is 
assumed to be small compared to the Andreev level spacing, V u -C E{t). 

The function u, which is smooth on the scale of the Fermi wavelength, has a discontinuity at x = 0, 
i.e. at the point contact, whose spatial extension can be neglected. The discontinuity translates to 
a boundary condition for u, 

u(+0) = (1/Vd)(1-VRt v )u(-0), R=l-D, (4) 

which can be found || by matching at the point x = two scattering state solutions to the BdG 
equation approaching from left and right. 

2.2. Resonant coupling of Andreev levels 

Under the assumption that eV, V g <C A the system of Andreev levels experiences an adiabatic 
evolution in time. This is true at all times except close to the resonances (points A and B in Fig. 1) 
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and at point C, which will be discussed later. For a resonant transition to occur, the deviation of the 
intcrlevel spacing from the resonance value E(t) — Tiuj/2 = E(tA,B)St, where St is the duration of the 
resonance, has to be less than the quantum mechanical resolution of the energy levels Ti/ St. Using this 
criterion we can estimate the duration St of the non-adiabatic evolution as St ~ [H/E^a^b)} 1 ^ 2 > which 
implies that St is much shorter than the period T p = frir/eV of Josephson oscillations if eV <C A. 
Hence, if this inequality is obeyed we may consider the non-adiabatic dynamics as temporally 
localized scattering events. We shall use this later to derive an analytical result for the induced dc 
current. For the time being, however, we will keep the discussion a little more general. We therefore 
introduce an Ansatz for the wave function u in the normal region of the junction in terms of a linear 
combination of the eigenstates corresponding to the adiabatic Andreev levels E ± . The rapidly 
oscillating terms are explicitly introduced in this so called resonance approximation, 

u(i) = b+ {t)e- wt / 2 u+ + b- {t)e lu)t/2 u- . (5) 

Following Ref H we insert the Ansatz (||) into the BdG equation (|J) to find, using the notation 
b = (b + ,b~), an equation for the coefficients 6 ± , 

vl_/n -5Lj+-{t) 

Here Sui-i (t) = [2E + (t) — hu>\/2Ti is a measure of the deviation from resonance. This equation 

describes the time evolution of the coefficients b + and b~ , which embody the dynamics of the pop- 
ulation of the Andreev states under irradiation. We recall that the coupling to the time-dependent 

field is finite only in the normal region, which explains why in (^|) obtains as the matrix 

element of a z V g (x,t) between the Andreev states u ± ; V-\ = (V u /2) J dx{u + cr z u~), where (..) de- 
notes a scalar product in 4-dimcnsional space. For the case of a double barrier SQPC structure 

V-\ was calculated in Ref. ||. For a single barrier junction an analogous calculation gives us, 

V-\ = a(L/£o) y/DRVu sin(</>/2), where the constant a ~ 1 is determined by the position of the 

barrier. We note, that this matrix element is proportional to the reflectivity of the junction; re- 
flection mixes electron states with +fc^ and —hp allowing optical transitions between the Andreev 
levels. In a perfectly transparent SQPC (D — 1), the upper and lower Andreev levels correspond 
to opposite electron momenta — the two levels cannot be coupled by radiation since momentum 
cannot be conserved — and the effect under consideration does not exist, cf. Refs. || (Tof . 

2.3. Boundary conditions at 4> = 27m 



ib(t) 



b(t). (6) 



Before we proceed to calculate the current through the irradiated SQPC we need to discuss the 
boundary condition at cf> = 2nn (n is an integer). In the vicinity of these points, the Andreev levels 
approach the continuum and the adiabatic approximation is unsatisfactory, even at small applied 
voltages and weak electromagnetic fields. The duration St of the non-adiabatic interaction between 
the Andreev level and the continuum states can be estimated using the same argument as for the 
microwave-induced Landau-Zener scattering. One finds that, St ~ h/(Ae 2 V 2 ) 1 ^ 3 . If we evaluate the 
condition, St << T p , we find that we would require that eV/A <C tt 3 . This means that we can 
safely treat the non-adiabatic region as short. To derive the boundary condition, for example, at 
point C in Fig. 1, one needs to calculate the transition amplitude connecting the states u + (ti) at 
time ti < t c - St and u+(t 2 ) at time t 2 > t c + St: (u + (t 2 ) U(i 2 , ii)u + (£i)). Here U(t 2 ,ii) is the 
exact propagator corresponding to the Hamiltonian in Eq. (j^) . We will now proceed with symmetry 
arguments to show that this amplitude is exactly zero. 

It can be shown that both the Hamiltonian (||) and the boundary condition for u at x = 
are invariant under the simultaneous charge- and parity inversion described by the unitary 
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operator A = P<j x t z , where P is the parity operator in x-space. This implies that at any time 
any non-degenerate eigenstate of the Hamiltonian is an eigenstate of the symmetry operator A 
with eigenvalue +1 or —1 and that this property persists during the time evolution of the state. 
Specifically, if we take a state on each side of <f> = 0, 

Au+(^) = A lU +(0) (7) 
Au+(-0) = A 2 u+(-<^). (8) 

Next we insert u + (—(f)) — Tu + (</>), where the operator T = a x r y changes the sign of the phase (f>, 
into Eq. (g) and apply a x r y from the left and arrive at, 

Au+(</>) = -A 2 u+(0), (9) 

which shows that Ai = — A 2 . This means that the two states, u + (0) and u + (—(f)) are orthogonal. 
This is consistent with the results of Shumeiko et.al. who have shown that the Andreev state 
wave functions are 47r-periodic whereas the energy levels and the current are 2-7r-periodic. Since 
the state evolving from the adiabatic state u + (ii) is orthogonal to the adiabatic state u + (i 2 ), the 
probability for an adiabatic Andreev state to be "scattered" into a localized state after passing the 
non-adiabatic region is identically zero. In reality, the Andreev state as it approaches the continuum 
band edge decays into the states of the continuum. Such a decay corresponds to a derealization in 
real space and is the mechanism for transferring energy to the reservoir [[fl| . 

The orthogonality property shown above guarantees that the coherent evolution of our system 
persists during only one period of the Josephson oscillation and that the equilibrium population of 
the Andreev levels is reset at each point </> = 27m |l2] ] . This imposes the boundary condition 

6+(27rn + 0) = 0, 6-(27m + 0) = 1, (10) 

at the beginning of each period. 

2.4. Calculation of the current 

The quasiclassical equation for the total time dependent current at the junction (x = 0) reads 
I(t) = —evp{\iT z u) \ X =Q. By manipulating the BdG equation as given by (||) and (|^) and omitting 
oscillating terms which will not contribute to the dc current — the quantity of interest — one may 
turn this expression into the following form, 

In the static limit, V g = and <f> — > 0, this result is clearly equivalent to the standard expression 
/ = (2e/K) (dE^/dcf)) for the Andreev level current. 

In the general nonstationary case u is a linear combination of u ± according to Eq. (|) and we 
can calculate the current averaged over the period T p . This is done through Eq. ( pi] ) with help of 
Eq.( ^J) and the normalization condition \b + \ 2 + \b~\ 2 = 1, leading to 

^ = ^{^-^)[\b + (T p )\ 2 -\bt\'], (12) 

where |& + (T p )| 2 is the population of the upper level at the end of the period and \b$\ 2 at the 
beginning of the period. The result simplifies further since according to Eq. ( |l(i| ) \b$\ 2 = 0. 

The direct current through the contact can be viewed as resulting from photon-assisted pair 
tunneling or equivalently as being due to the distortion of the ac pair current due to the induced 



d(uii) • 
ih— V g (ua z u) 



(11) 
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intcrlcvel transitions. To facilitate an understanding of the current expression in another way, let us 
study energy conservation. We have two sources of energy, the applied field and the applied bias. If 
we consider a single Josephson period T p , the energy absorbed by the system from the voltage source 
is I ■ V ■ T p , and from the applied field it is fiuo\b + (T p )\ 2 (assuming for simplicity that |fyj~| 2 = 0). 
Energy conservation for the period can be stated as, E out = E a ^ sor b e ^, with the energy leaving the 
system stated as 2A|6+(T p )| 2 : 

2A|6+(T P )| 2 = /• V ■ T p + htu\b + {T p )\ 2 . (13) 

The current can then easily be found as, 

I = J_(2A|6+(T p )| 2 - M& + (T p )| 2 ) (14) 

or, inserting the value of T p , 

/ = £ (A ~T )|6+(Tp)|2 - (15) 

This is exactly the same result as obtained by the more tedious method used above. This discussion 
allows us to interpret the current as voltage-bias mediated, photon assisted tunneling. 



3. Results of numerical and analytical analysis 

Generally the current is given by Eq. ( p"2| ) with the boundary condition Eq. ([To]) inserted. The actual 
calculation of the current has to be done numerically in most cases. However, in the limit of a weak 
applied external field and when the frequency of the applied field is such that we have two (in time) 
well separated resonances, we can treat the resonances as temporally localized scattering events. 
This approach can be quite rewarding as will become clear below. 



3.1. Scattering formalism 



Formally, we can describe the system's evolution through a resonance by letting a scattering matrix 
S connect the coefficients b before and after the splitting points A and B. Approximating the time 
dependence of the Andreev levels to be linear, a standard analysis of the Landau-Zener interlevel 
transitions (see e.g. ra) gives the scattering matrix elements at the points A and B as, 



S A 



r 

de ie 



-de 



-is 



, Sb = 



r 

-de- 10 



de 



(16) 



where |<i| 2 



1 



1 — e 7 is the probability of the Landau-Zener interlevel transition. Here 



7 = tt\V-\ \ 2 /\dE/dt\ where Vj is once again the matrix element for the interlevel transitions. 

By introducing the matrix = exp (ia z $(i, j)), 



1 



2eV 



d0 E(cf>) 



hto\ 



(17) 



which describes the "ballistic" dynamics of the system between the Landau-Zener scattering events, 
we connect the coefficients 6 ± at the end of the period of the Josephson oscillation, 4> — 0, with the 
coefficients b^ at the beginning of the period, 4> = —2-k, 



&0,bSb&B,aSa$A,-2tt I , _ 



(18) 




Fig. 2. Current vs. inverse voltage from Eq. (15) for a biased SQPC irradiated with microwaves 

of frequency oj = 1.52A/7L and amplitude corresponding to a matrix element \V-\ | = 0.024A for 

interlevel transitions, I c = eAD/ti. Note the cut in the inverse voltage scale. Results of the scattering 
approach, Eq. (19) [o], are close to those obtained by numerically solving the BdG-equation with 
the radiation field treated in the resonance approximation (solid line). The close fit means that 
the scattering picture can be used to reconstruct the Andreev levels from the period of the current 
oscillations (Andreev level spectroscopy, see Section 3.2). 



The time-averaged current through the junction can be directly expressed through these coefficients. 
Combining the boundary condition (|l^) with Eqs. ( |l2| ) and (|l|), one finds 

I dc = ^(a-^-) e -T(l - e -f) sin 2 ($(A, B) + 9), (19) 



ttH \ 2 r 

where 9 is the phase of the probability amplitude for the Landau-Zener transition, which only weakly 
depends on V . In Fig. 2 the current is calculated both numerically and with the expression above, 



Eq. (19). The current is plotted as a function of inverse voltage in order to show a periodicity in 
1/V. The close fit in Fig. 2 ensures us that the scattering approach works well for the situation 
when the resonances are well separated and the applied field is weak. 

3.2. The quantum point contact as an interferometer 

Equation ( |l9| ) is the basis for presenting the biased SQPC as a quantum interferometer. There is 
a clear analogy between the SQPC interferometer and a standard SQUID in that they both rely 
on the presence of trajectories that form a closed loop. In a SQUID, which is used to measure 
magnetic fields, the loop is determined by the device geometry; in the SQPC the voltage (analog of 
the magnetic field) is well defined while the "geometry" of the loop in (E, t) space can be measured. 
This loop is determined by the Andreev-level trajectories in (E, t) space and is controlled by the 
frequency of the external field. This gives us an immediate possibility to reconstruct the phase 
dependence of the Andreev levels from the frequency dependence of the period II of oscillations of 
the current versus inverse voltage, see Fig. 2. Indeed, it follows from Eq. (0), that 



dh> 



(20) 

IE 



This shows that the position of the Andreev energy levels can be reconstructed from measurements 
of microwave-induced subgap currents through the SQPC. 

4. Discussion - dephasing and relaxation 



In order to be able to do interferometry it is necessary to keep phase coherence during at least one 
period of the Josephson oscillation. There are three dephasing mechanisms that impose limitations 
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in practice: (i) deviations from an ideal voltage bias, (ii) microscopic interactions, and (iii) radiation 
induced transitions to continuum states. The main source of fluctuations of the applied voltage is 
the ac Josephson effect. According to the RSJ model, a fixed voltage across the junction can only 
be maintained if the ratio between the intrinsic resistance Ri of the voltage source and the normal 
junction resistance Rn is small. If Ri/Rn <§; 1 the amplitude of the voltage fluctuations 5V is 
estimated as SV ~ (Ri/ R^)A. Effects of voltage fluctuations on the accumulated phase &(A, B) 
can be neglected if 6$(A, B) = &(A, B)SV <C 2n, i.e. if eV > (Ri/ 'Rn) 1 ^ 2 A, which corresponds to 
a lower limit for the bias voltage [Q . 

The dephasing time due to microscopic interactions is comparable to the corresponding relaxation 
time This mechanism of dephasing can be neglected as soon as the relaxation time exceeds 
the Josephson oscillation period, Tj 3> hn/eV. Taking electron-phonon interaction as the leading 
mechanism of inelastic relaxation, we estimate Tj to be of the order of the electron-phonon mean 
free time at the critical temperature, T p h(T c ), since the large deviations from equilibrium in our case 
occur in the energy interval E < A. This gives fl5|| another limitation on how small the applied 
voltage can be, eV > 10 _2 A. 

The third mechanism of dephasing becomes important when the Andreev levels are closer than 
huj to the continuum band edge. One can estimate the corresponding relaxation time as r w ~ 
hA/V 2 . For small radiation amplitudes r w exceeds Ty, while for optimal amplitudes they are about 
equal. The effect of the level-continuum transitions on the interference oscillations depends on 
the frequency. If too < 2A/3 the "loop region" [E(t) < hu>/2] is optically disconnected from the 
continuum, and transitions cannot destroy interference. Possible level-continuum transitions at times 
outside the loop will only decrease the amplitude of the effect by a factor exp(— aV 2 /eVA), where 
a < 1 is the relative fraction of the period Ty during which transitions to the continuum are possible. 
Accordingly, this factor is of the order of unity for the voltages that correspond to the maximum 
amplitude of oscillations. If Hlu > 2 A/3, the interference is impeded by the optical transitions into 
the continuum, and the current oscillations decrease. Still, a nonzero average current through the 
junction will persist. 

In any real system, both relaxation and dephasing will be present as discussed above. Therefore, to 
complete our analysis we will simulate these effects on our system. If we assume that the relaxation 
and dephasing times, t,t^ are long compared to the duration of the non-adiabatic resonances, 5t, 
we can use a technique JL4|, where dissipation is modelled by adding a term to the time evolution 
equation of the density matrix for the two-level Andreev system. The equations are, where relaxation 
enters in the diagonal terms and dephasing through the off-diagonal terms, 

Pnn(t) = - l -[H Q (t),p(t)] nn - (21) 

Pnn>(t) = -U Ho (t),p(t)) nn , - n^n'. (22) 

n T<f, 

with r as the characteristic time for relaxation of the system to the equilibrium population, p eq and 
Ttj) as the characteristic time for dephasing. These equations describe the "ballistic" dynamics of the 
dissipative system, replacing Eq. (|l7|). 

The exact form of the density matrix for the system considered here will be a 2x2 matrix for 
the discrete two level space of the Andreev states. (To avoid confusion we will use <7i,<72,<73 to 
denote the Pauli matrices in this discrete space.) The Hamiltonian for the density matrix is found in 
Eq. ([|), and the diagonal elements p\\ and pii will represent the population of the upper and lower 
Andreev levels. The resonant events which are unaffected by the dissipation can still be modeled 
with the scattering matrices introduced earlier, Eq. (|l6|), 

p' = S A ,BpS\ B . (23) 
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To calculate the effect of relaxation and dcphasing on the current we will use the standard 
expression, I — 2e/h(dE/d(f)) ) for the current carried by a populated Andreev state and average 
over one Josephson period. The result is, as a function of the upper levels population, 

= Y p jf ' = T p tC T ^3p)9£/50 dt = -L £" [ 2pu - 1] dt, (24) 

where we have used the Josephson relation <fi = 2eV/h and the conservation of probability, Tr{p) = 1. 

By solving the time evolution equations of the density matrix for the whole Josephson period and 
imposing the boundary condition p(Q) = p eq = (1 — cr^)/2 ,we arrive at the following result, 



Pn = < d 2 



,-(*-ti)A 



, < t < ti 
,t 1 <t<t 2 



/(r 2 ,r ,r) e"^-* 2 )^ , t 2 < t < T p 



where 

/(r 2 ,r ,r) 



d\l 



-(ta-ti)/n 



2r 2 d 2 (e 



-(t 2 -ii)/r 



-(ta-ii)/T> 



(25) 



cos2($(A,B) + e)) . (26) 



Inserting these expressions into Eq. (p3) we find, 



jdiss 



(V) 



d 2 hu> 



^ e -(ta-ti)/r _ + f 2 E +(t)dt+ 



/(r 2 ,r ,r) 



\ e -(T p -t 2 )/r _ ^ 

2 



E+(t)e- (t - t2)/T dt 



and 



it reduces to the 



The first important property of this expression is that for - 
previous current expression that we derived, Eq. (|l^). One obvious effect on the current is that the 
oscillations of the current decrease when the dephasing time, r^, becomes small. Another is that 
when the relaxation time, r, is decreased, the current diminishes, see Fig. 3. 

We note that the number of parameters has become quite large. It is now possible to calculate 
the current for different combinations of r, r^, V, V-\ , us, and the transmission coefficient D. 

A curious effect occurs when the relaxation time is such that the upper level, after being populated 
at point A, is fully depopulated at = 7r, the middle of the Josephson period. What we find is a 
negative dc current for a positive bias, a regime where our system shows a negative conductance, 
see Fig. 3, plots a and d. The relaxation time is such that the upper level is mainly populated when 
it's derivative is negative (it carries a negative current) and then the lower level is highly populated 
for the part of the period when it's derivative is negative. This effect will be most pronounced when 
the energy of the applied field, Huj, approaches the energy of the gap, 2A, as in Fig. 3d. 

Actually, the physical mechanism which is behind the appearance of this negative resistance is 
similar to the one responsible for the "somersault effect" effect discussed in Ref. ||. In that case 
the method for periodically bringing a pair of Andreev levels in a superconducting quantum point 
contact in resonance with a microwave field was different. Rather than applying a voltage bias, the 
time-dependence of the gate potential was assumed to have two components; one slow component 
which shifts the Andreev levels through its effect on the transparency D — D(kp(t)) and one fast 
(microwave) component for coupling the two levels. 



5. Conclusions 



We have shown that irradiation of a voltage-biased superconducting quantum point contact at 
frequencies u> ~ A can remove the suppression of subgap dc transport through Andreev levels. 
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Fig. 3. Simulations of the current when relaxation and dcphasing are added. In a) only relaxation 
is added, with r sa 2 • 10~ 10 s, the average value of the current decreases and also may become 
negative. In b) there is only dephasing, with w 5 ■ 10~ 10 s. In c) both types of dissipation are 
added with r = w 5 • 10~ 10 s. When, as in d), the frequency of the applied field is increased such 
that it approaches the superconductor energy gap, 2A, we find that the system shows a negative 
conductance, see text, when r = sa 2 • 10 _10 s. For plots a,b,c hco/2A w 0.75 and for plot d 

huj/2A w 0.95. The amplitude of the external field is such that | |/A = 0.024 in all plots, and 

for reference the case when r — * oo and — •> oo are included (dotted curves). 



Quantum interference among resonant scattering events can be used for microwave spectroscopy of 
the Andreev levels. 

The same interference effect can also be applied for detecting weak electromagnetic signals up 
to the gap frequency. Due to the resonant character of the phenomenon, the current response 
is proportional to the ratio between the amplitude of the applied field and the applied voltage, 
I ~ \V±\ /AeV. At the same time, for common SIS detectors a non-resonant current response is 
proportional to the ratio between the amplitude and the frequency of the applied radiation 
I ~ |V±/c!j| 2 , i.e. it depends entirely on the parameters of the external signal and cannot be improved. 

Finally, we note that the classic double-slit interference experiment, where two spatially separated 
trajectories combine to form an interference pattern, clearly demonstrates the wave-like nature of 
electron propagation. For a 0-dimensional system, with no spatial structure, we have shown that a 
completely analogous interference phenomenon may occur between two distinct trajectories in the 
temporal evolution of a quantum system. 
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